Universality, the QCD critical/tricritical point and the quark number susceptibility 
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The quark number susceptibility near the QCD critical end-point (CEP), the tricritical point 
(TCP) and the 0(4) critical line at finite temperature and quark chemical potential is investigated. 
, Based on the universality argument and numerical model calculations we propose a possibility 

£SJ ' that the hidden tricritical point strongly affects the critical phenomena around the critical end- 

— » | point. We made a semi-quantitative study of the quark number susceptibility near CEP/TCP for 

several quark masses on the basis of the Cornwall-Jackiw-Tomboulis (CJT) potential for QCD in 
the improved-ladder approximation. The results show that the susceptibility is enhanced in a wide 
region around CEP inside which the critical exponent gradually changes from that of CEP to that 
of TCP, indicating a crossover of different universality classes. 

(N ' 
> 

The vacuum of the quantum chromodynamics (QCD) is believed to undergo a phase transition to the quark-gluon 
plasma (QGP) at high temperature T and/or at high quark chemical potential [i. Such a new state of matter is 
expected to be produced in on-going heavy-ion collision experiments at Relativistic Heavy-Ion Collider (RHIC) and 
^vq , in the future Large Hadron Collider (LHC) [1]. 

The phase transition of the hadronic matter to QGP at finite T with /i = has been studied extensively on the 
lattice. In particular, the chiral phase transition is likely to be of second order for QCD with two massless quarks. Also, 
the static critical behavior is expected to fall into the universality class of the 0(4) spin model in three dimensions [2]. 
In nature, the light quarks have small but finite masses and the second order transition becomes a smooth crossover. 
Oh, Study of the QCD phase transition with finite /i has been retarded because reliable lattice simulations have not 
been available so far due to the severe fermion sign problem. Nevertheless, there is a growing evidence that the phase 
diagram of QCD with massless 2-flavors has a tricritical point (TCP, Fig. 1, point P) at which a line of critical points 
(the 0(4) line) at lower /i's turns into a first order phase transition line at higher /i's. The existence of TCP was in 
fact suggested in various calculations based on effective theories of QCD [3-10]. If the u, d-quark masses are increased 
from zero, a line of critical points (the wing critical line) emerges from TCP and the point which corresponds to the 
physical quark mass m p hy S is called the QCD critical end-point (CEP, Fig. 1, point E) because this is the point where 
the first order phase transition line terminates. Indeed, some evidence of the existence of CEP was shown recently in 
a lattice QCD simulation with 2+1 flavors by Fodor and Katz [11]. In this paper we assume that CEP exists in the 
phase diagram of QCD. 1 

Second order phase transitions are characterized by the long- wavelength fluctuations of the order parameter. In 
the case of CEP, it is the sigma (a) field. Then, it is expected that the fluctuations of the sigma field will be reflected 
in the event-by-event fluctuation of pion (ir) observables due to the the ir — a coupling. Based on this observation, 
possible observable signals associated with CEP have been studied in detail in relation to the relativistic heavy-ion 
collision experiments [12-14]. 

The purpose of this paper is to point out that the anomaly near CEP is not pointlike but has much richer structure. 
Our starting point is a simple question: "How large is the critical region?" The critical region is defined as the region 
where the mean field theory (or the Landau theory) of phase transitions breaks down and the true non-trivial critical 
exponents can be seen. Usually, one expects that the critical region is surrounded by the mean field region and the 
critical exponents change from the non-trivial values to the mean field values as one comes away from the critical 
point. One might argue that this question is only of academic interest because the non-trivial exponents and the mean 
field exponents are numerically not so different and probably experiments cannot distinguish them. [This observation 
is the basis of [12].] However, as we will see, pursuing this question leads to an important notion which may shed light 
on certain results of both heavy-ion collision experiments and future lattice simulations at finite chemical potentials. 



1 Accordingly, we fix the strange quark mass to its physical value. Below 'the quark mass' means the u, d current quark mass 
which we consider as a variable parameter. 



There is a well-known criterion which estimates the size of the critical region, the Ginzburg criterion [15]. It tells 
that if the singular part of the thermodynamic potential f2 (the Landau-Ginzburg potential) for a certain second order 
phase transition is given by 



^singular = c(V (j)) 1 + Clt(f) 2 + b(f> 4 , (1.1) 

where <j) is the order parameter and t = T ~ Tc is the reduced temperature (T c is the critical temperature in the mean 
field approximation), the critical region is estimated to be 

T 2 b 2 

1*1 ~^r- d-2) 

At first sight, this criterion seems useless because we do not know the coefficients appearing in (1.1) for CEP. 2 However, 
in the next section we will derive a bound to the size of the critical region. In fact, there is a reason to expect that 
the critical region of CEP is small. This is because the QCD critical end-point is a descendant of the tricritical point 
of the massless theory. 

This observation led us to study the critical phenomena of both CEP and TCP simultaneously and their possible 
correlations. We make both qualitative and quantitative analyses of the physics near TCP and CEP with particular 
emphasis on the (singular) behavior of the quark number susceptibility Xq defined by 

ld 2 VL , , 

x « = -vW (L3) 

where 0(T, f£) is the thermodynamic potential and V is the volume of the system. Xq lSi a response of the quark 
number density to the variation of the quark chemical potential and is one of the key quantities characterizing the 
phase change from the hadronic matter to QGP [18-22]. The lattice data tell us that, at fi — 0, Xq increases rapidly 
but smoothly near the critical temperature [23-25]. On the other hand, the universality argument predicts that it 
diverges at both TCP and CEP with certain critical exponents. Therefore, it would be important to study its critical 
behavior with and without the quark masses to see whether it can provide a new way of detecting the TCP / CEP on 
the lattice as well as in the heavy-ion collision experiments. 

In addition to Xq, we occasionally mention the singular behavior of the specific heat C and the chiral susceptibility 
Xch defined as 

c = -vdf 2 > (L4) 
i d 2 n , , 

u= w (L5) 

From the viewpoint of critical phenomena, Xq an d C are essentially the same near TCP/CEP while Xch is different 
from Xq near TCP and only slightly different near CEP in the sense that will be clarified below. 

In Section II, we will make a general analysis of the interplay between TCP and CEP in the small quark mass limit 
based on the universality argument. After determining the relative location of TCP and CEP in the phase diagram 
as a function of the quark mass, we construct the Landau-Ginzburg potential for CEP to determine the singular 
behavior of susceptibilities both in and beyond the mean field approximation. It turns out that the smallness of the 
quark mass gives a bound to the growth of the critical region, turning our attention to the tricritical point. Then we 
discuss a possible crossover from the tricritical universality class to the Ising universality class. 

The universality argument is so general that it gives no quantitative results. In order to reinforce the ideas given in 
Section II, in Section III we show the results of the numerical calculation on a model, the Cornwall- Jackiw-Tomboulis 
(CJT) potential for QCD [26] in the improved ladder approximation [10,27]. We will find how well the numerical 
results match with the qualitative predictions of the universality argument, demonstrating the power of universality. 



2 The size of the critical region depends on the microscopic dynamics and universality tells nothing about it. A clear example 
is the A-transition of liquid helium and the superconducting transition of metals. Although they belong to the same universality 
class (the 0(2) spin model), their critical regions are very different; \t\ ~ 0.3 for the A-point and \t\ ~ 1CT 15 for conventional 
(type-I) superconductors. Just for reference, we note that for typical liquid-gas phase transitions which belong to the same 
universality class as the phase transition at CEP, \t\ ~ 1CP 2 [16]. [Corrections to the scaling [17] are not negligible until one 
reaches \t\ ~ 1(T 4 .] 



In particular, we observed some indication of the effect of TCP on the QCD phase diagram even with a reasonable 
value of the quark mass. 

Section IV is devoted to conclusions. 

Brief description of the model is given in Appendix A. In Appendix B, we discuss, for completeness, how Xq behaves 
along the 0(4) line based on the universality argument. We will see that the monotonous increase so far observed on 
the lattice is a property only at \i = 0. 

II. UNIVERSALITY ARGUMENTS 

Universality is such a strong notion of modern physics [28] that its applicability ranges from phase transitions in 
ordinary liquids to thermal phase transitions of relativistic quantum field theories. In this section we study the critical 
phenomena near CEP/TCP based on the universality argument. We will see that a lot of general information can be 
extracted by the universality argument alone without mentioning any complexities of the strong interaction. 

A. The QCD critical end-point 

It was suggested theoretically [4,5,8,10] and found on the lattice [11] that QCD has the critical end-point (CEP) 
at finite temperature T c and baryon chemical potential ji c (Fig. 1, point E). At the critical end-point, only the 
cr-field becomes massless and the universality class of this phase transition is considered to be the same as that of the 
liquid-gas phase transition, or equivalcntly, that of the 3-dimensional Ising model. 3 

In order to exploit the power of universality to investigate the singular behavior of various quantities, we consider 
the mapping of the (ti,hi) axes of the Ising model (ti is the reduced temperature and hi is the reduced magnetic 
field) onto the (T,n,m) space (m is the light quark mass divided by the typical scale of the problem such as T c ). 
This can be achieved by considering the tricritical point (TCP) at (T,^,m) = (T t ,^ t ,0) (Fig. 1, point P). Below we 
explicitly construct the Landau- Ginzburg potential for CEP starting from the general theory of tricritical points [29] 
and discuss associated universal behaviors. 



3 This is not obvious a priori and requires explanation. As we shall see below, the phase transition at the end-point is 
characterized by the one-component order parameter a. The effective Landau- Ginzburg potential contains odd powers of 
a which break the a — > —a symmetry of the Ising model. This is the same situation as the liquid-gas phase transition. 
Theoretically, the usual renormalization group argument should be reconsidered in the presence of the asymmetry [30] . Although 
there are some subtleties about this problem, experimentally it is clear that the liquid-gas phase transition and the 3D Ising 
model belong to the same universality class. 
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FIG. 1. The phase diagram of QCD in the (T, fi, m) space. Point P is the tricritical point (TCP) of the massless theory and 
point E is the critical end-point (CEP) of the real world. The dotted lines represent the second order phase transition and the 
solid line represents the first order phase transition. 



Near TCP, long-wavelength physics of the system can be described by the thermodynamic potential expanded up 
to the sixth-order in the order parameter field (the sigma field) a 

n M F = ^o-mcr+^a 2 + ^a' l + ^a 6 , (2.1) 
2 4 o 

where flo is the contribution from short wavelength degrees of freedom irrelevant to the study of critical phenomena. 

At TCP, a = b = m = 0. Assuming that a and b are analytic in T and fi and that c > is approximately constant 
near the tricritical point, we expand them as follows [31] 

a(T,ii)=C a (T-T t ) + D a (ji-iH) 

b{T,n)=C b {T-T t ) + D b (n-n t ), (2.2) 

where we have neglected terms higher order in the deviation from the tricritical point. C a > and D a > are related 
such that the line a(T, n) = is tangential to the first order phase transition curve at TCP. b is positive for T — T t > 
(fj, < m) on the a = line, which leads to the condition 

C b D a - C a D b > 0. (2.3) 

These conditions come from the geometry of the phase diagram, namely, the fact that there is a line of (bi-)critical 
points at T > T t , \l < [i t - We do not know the actual values of these coefficients. But we need not know them for the 
present purpose. 

If we increase m from zero, at some point (T c (m), /z c (m)) in the (T,fi) plane two minima and a maximum of the 
potential coalesce. This is the critical end-point. There the sigma field acquires a non-zero expectation value cr a which 
is determined by the following equations [In this section we exclusively consider the small m limit and leave only the 
leading terms in to.] 

n'(T c (m),/i c (m), <7 ) = -to + a m <7 + b m al + cerfj = 
fl"(T c (m),fi c {m),a ) =a m + ib m ul + 5ca^ = 

Q'"(T c (to),McK), <r ) - 6b m a + 20 C( r 3 = 0, (2.4) 



where a m = a(T c (m) , /U c (to)) and b m = (T c (m) , fi c (mj) . The solution is 
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(2.5) 



Using (2.2) and (2.5) we can locate the critical end-point for small m. 
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(2.6) 



(54)*(C 6 D -C D 6 ) 

Thus, as we increase the quark mass m, the critical temperature decreases and the critical chemical potential increases 
at least for small m [Fig. 1]. Expanding fi(T, fi, a) around (T c (m) , n c (m) , ao) we obtain the Landau-Ginzburg potential 
with the new order parameter a = a — gq 

&mf{T,Ii,o) = tt MF (T c (m),fj, c (m),(T ) 

+ A 1 a + A 2 a 2 + A 3 a 3 + A i a 4 , (2.7) 



where 



At 



A 2 



A 4 



(C a a + C b o%)(T-T c (m)) 
+ {D a ao + D b crl)(n - (J, c (m)), 
(C a + 3C b a 2 )(T-T c (m)) 
+ (D a + 3D b o$)(ii- fj, c (m)), 
C b (T - T c {m)) + D b (ji - n c {m)), 
-b m 



(2.8) 



Ai (i = 1,2,3) vanish at the critical point whereas A4 does not, indicating that (T c (m) , ^i c (m)) is an ordinary 
(bi-)critical point as stated above. 

Looking at (2.7) and (2.8), we immediately notice two important things. First, A 2 is not proportional to T — T c (m) 
as in the Landau theory but a linear combination of T — T c (m) and /1 — /i c (m). This means that T and fj, are equivalent 
thermodynamic variables in the sense of Griffiths and Wheeler [32] and that A 2 is the temperature-like scaling field 
which corresponds to ti of the Ising model. Second, A\ rather than the quark mass plays the role of the 'external 
field' which is conjugate to the new order parameter. Thus it can be identified as the magnetic field-like scaling field 
hi. Indeed, it is easy to show that, on the line A\ — 0, A 2 and A 3 are positive for T > T c (m) (or \i < fj, c (m)) and 
negative for T < T c (m) (or /x < fi c {m)) and this line is asymptotically parallel to the first order phase transition line 
at the critical end-point. See Fig. 2. 




FIG. 2. The mapping of the Ising model axes onto the (T, fx) plane. The solid lines represent the first order phase transition 
(the coexisting line). The dashed lines separate regions with different exponents. 



Now we can discuss the critical behavior of susceptibilities; the quark number susceptibility \ q , the specific heat C 
and the chiral susceptibility Xch- I n the mean field approximation, the equilibrium value of a is determined by the 
first and fourth order terms of (2.7) in the small mass limit. Then we obtain, for paths asymptotically not parallel to 
the A\ = line (the first order phase transition line) 

Xq ~ C ~ mfs\g - ,g c |~ e , 

Xch ~m-T§ \g-g c \~ e , (2.9) 

where e = = |. \g — g c \ denotes the distance from CEP in some units. For the path asymptotically parallel to the 
Ai = line, the exponent is 7 = 1 > e. Note that, although the critical exponents are the same, the amplitude of the 
chiral susceptibility is enhanced whereas that of the quark number susceptibility is suppressed by factors of m. 

Inside the critical region, where the mean field theory breaks down, Q\ singular does not admit a simple expansion 
with smooth coefficients. (2.7) should be regarded as the saddle point approximation to the following functional 
integral 

Sl(T,tM,a) = -fin I [da , ]exp(~ J d 3 rH eff (r)), (2.10) 

where H e ff is the Landau-Ginzburg- Wilson hamiltonian 

H eff = A' (Vcr r ) 2 + A\a' + A' 2 a 12 + A^cr 13 + A' 4 a' 4 . (2.11) 

A\ (i = 1 ~ 4) arc in general different from Ai due to fluctuations. However, we expect that the differences between 
A\ and Ai are of the higher order in m. 4 Note the appearance of the kinetic term. The sigma field is no longer a 
constant beyond the mean field approximation. The potential (2.10) will eventually lead to the scaling equation of 
state [33] written in terms of the scaling fields A\ and A2 (the revised scaling [34]). Because T, /1 and m participate 
in the magnetic field-like scaling field, we obtain, very schematically, the most singular part 5 

d 2 n 2 d 2 n 2 . . .. 

X ^W~ m ~ 5 dAj~ mHaa) 

d 2 n 2 _d 2 n 2 , 
d 2 n _ ± d 2 n 4 . . 

Xch ~dnT^ m ^dAj ~™~^> 

(aV) ~ \g-g c \- £ (2.12) 

where \g — g c \ is the distance from the critical end-point in some units, e = 7/ [35 w 0.8 for any direction which make 
an angle to the A\ = line at the critical end-point. For the path asymptotically parallel to that line, the exponent 
is 7 <~ 1.2 > e. [These values are taken from the 3D Ising model.] 

Having discussed the singular behavior of susceptibilities inside the critical region, however, we give a pessimistic 
result. Since we now have the Landau-Ginzburg potential for CEP, we can say something about the size of the critical 
region. Recall that, according to the Ginzburg criterion (1.2), the radius of the critical region is proportional to the 
square of the coefficient of the quartic term. Other coefficients are quark mass independent in the leading order. Thus 
we obtain 

|i|~A|~mt. (2.13) 

This gives a bound to the size of the critical region. It shrinks to zero as the quark mass decreases (See, Fig. 1). The 
physical reason behind this is that the coefficient of the quartic term is zero at the tricritical point and remains small 
near it. 



4 The coefficients are further affected by the change of integration variables. These degrees of freedom can eliminate A' 3 , but 
do not change A' 12 in the leading order. In fact, only the direction of A' x is important to discuss the behaviors of quantities 
considered here (i.e., second derivatives of Q in directions parallel to the T, fi and m axes) [32,34]. 

5 In calculating Xch, dominant contribution to ££± comes from <r a d< ~ T ~J^ m ^ rather than ^^(T — T c (m)). The latter term, 
being proportional to T — T c (m), behaves as the correction to the scaling. Also, if the derivative acts on A2, we get (<j' 2 a' 2 ) ~ 
\g — g c \~ a ^ 5 which is less singular than the magnetic susceptibility. 



Generally speaking, the critical point of a strongly interacting system has a large critical region [35] . Thus the size 
of the critical region of CEP is subject to a competition between these opposite effects and the determination of it is 
a highly nontrivial problem. However, it seems to us that the above bound (2.13) is a compelling reason to expect 
that the critical region is 'small'. 

If the critical region of CEP is small, probably most of the fluctuations associated with CEP come from the mean 
field region around the critical region. 6 The central point of this paper is that if we consider the mean field region 
belonging to CEP, we should also consider the mean field region belonging to TCP. The tricritical point has, so to 
speak, a 'tricritical region' (see Fig. 1) which is a sphere or an ellipsoid in the (T,/i,m) space centered at (T t ,/z t ,0). 7 
Then it is possible that the tricritical region survives in the physical (T, /i) plane. The magnitude of the u, <i-quark 
masses is crucial to this. More interesting possibility is that the critical point is inside the tricritical region and a 
crossover of different universality classes happens (not to be confused with the crossover phase transition at lower 
chemical potentials) . Namely, as we approach CEP the critical exponents gradually change from those of the tricritical 
point to those of the 3D Ising model via those of CEP in the mean field approximation. [Note that the mean field 
exponents of a bi-critical point are different from those of a tricritical point.] Indeed, such kind of crossover was 
experimentally observed in an antifcrromagnct dysprosium aluminium garnet (DA1G) long ago. The critical exponent 
(3 for the magnetization tends to change from the tricritical value (/? — 1) to the Ising model value (/3 ~ 0.31) as we 
go along the wing critical line [36]. 

Thus, through the consideration of the critical region, we have become aware of a possible interesting role played 
by the hidden tricritical point. Its critical phenomena are therefore worth studying and will be discussed in the next 
subsection. 



B. The QCD tricritical point 



Motivated by the above arguments, we now turn our attention back to the QCD tricritical point. Because the 
upper critical dimension of models described by (2.1) is three, the origin of the coupling constant is an attractive IR 
fixed point. Correspondingly, universal behaviors associated with the tricritical point are well described by the mean 
field theory up to logarithmic corrections. 8 

Let us see how susceptibilities scale with respect to T — T t |, \fi — /x t | and m in the mean field approximation. 

At (T, (i, m = 0), straightforward calculations show that 

Xq ~ \h-h t \-~<\ 

Xch~\h-h t \-~«\ (2.14) 

where \h — h t \ is the distance (in some units) from TCP in the (T,fi) plane. 7 9 = |, "f c h = 1 for paths which arc not 
asymptotically tangential to the first order phase transition line. 

At (T,|j,m/0), the expectation value of a is given by the following equation 

m = acr + ba 3 + a 5 . (2-15) 

Near (T t ,^t,rn) where a = b = 0, [Note that this is the 'nearest' point to TCP in the phase diagram with a quark 
mass m.] we can expand the solution up to the second order in a and b 

i a _3 b _i o _7 o _3 i. . . 

(j = m 5 m 5 m 5 + U(a m 5 ,bm 5 ,abm ). (2-16) 



It must be cautioned that the mean field region does not always exist. For example, it is known that there is no mean field 
region for the A-transition of liquid helium (the critical region is large, \t\ ~ 0.3). However, if the critical region is squeezed 
by an explicit parameter of the theory as in the present situation, it would be meaningful to discuss the mean field region 
belonging to the critical point. [We thank M. A. Stephanov for a discussion on this point.] 

7 Here we use the term 'tricritical region' loosely for the region where any mean field-like effects of the tricritical point on 
susceptibilities exist. This terminology is a bit misleading because there is no critical region for a tricritical point in the usual 
sense. 

s This is why we neglected the pion degrees of freedom in (2.1). Mean field theory is truly universal in the sense that it does 
not depend on even the symmetry of the order parameter. However, the multiplicative logarithmic corrections to the scaling 
do depend on the symmetry of the order parameter. 



Inserting (2.16) into (2.1) and differentiating with respect to fi twice, we get \q- Because of (2.2) the differentiation 
with respect to \i is replaced by the differentiation with respect to a and b. Extracting the most singular contribution, 
we obtain 

Xq T^-\a=b=0~m 

Analogously, 

Xch~m-%. (2.18) 

The divergence of Xq is rather moderate in the mass direction, from which we expect that the quark number sus- 
ceptibility may still be large even with non-zero quark masses. Indeed, from (2.9) and (2.13) we can derive the 
ra-dependence of Xq at the edge of the critical region 

2 2 2 4 2 2 

X q ~ m~j |i|~3 ~ mis (ms) _ 3 ~ m~~ . (2-19) 

Comparing with (2.17), we see that the m-dependence is exactly the same. There may or may not be a reason for 
this coincidence. In any case, this does show that TCP is as important as CEP at least in the small quark mass limit. 

Starting from the simple Landau- Ginzburg potential, we have extracted a lot of physics near CEP/TCP. These 
analyses show the power of universality as well as its limitations. For example, the universality argument does not 
tell us whether the effect of TCP survives in the (T, //) plane with the quark mass of, say, 5 MeV. In order to quantify 
the ideas given in this section, we must resort to a specific model. This will be the subject of the next section. 



(2.17) 



III. NUMERICAL RESULTS 



In this section, we numerically calculate the quark number susceptibility in the (T, n) plane by using a model. As 
expected, the susceptibility diverges both at the critical and tricritical points. We also calculate the corresponding 
critical exponent. The results clearly demonstrate that the hidden tricritical point can affect the phase diagram with 
nonzero quark masses. 



A. CJT effective potential and the chiral phase transition 



As a model, we employ the Cornwall- Jackiw-Tomboulis (CJT) effective potential [26] for 2-flavor QCD in the 
improved ladder approximation [10]. A brief description of the model is given in Appendix A. For more details, see 
[10]. 

At zero temperature and chemical potential, the effective potential V is given by 

dp 2 p 2 

where the gauge coupling constant g 2 (p 2 ) and the dynamical quark mass function E (p) are 

»V) = 1 u u 2l l 2T > (3-2) 

ln[{p 2 + pi) /Aqqd J 

S(p) = m q (ln[(p 2 +^)/A QCD 2 ])- a/2 + ^— (ln[(p 2 + pD/A^ 2 })^ 1 • (3-3) 

p ~r P c 

p c is a momentum scale which separates the infrared (nonpcrturbative) region from the ultraviolet (perturbative) 
region. Cp = {N 2 — l)/2N c is the quadratic Casimir operator for the fundamental representation of the color SU(iV c ) 
group and a = 24/ (lliV c — 2Nf) (Nf is the number of active flavors. 9 a is proportional to the renormalization group 



'Although the potential is evaluated with Nf — 2, we take Nf — 3 in the gauge coupling (3.2). In this way we include the 
effect of s-quark only through the running of the coupling constant. 



invariant chiral condensate (qq) as a = 2n 2 a(qq) /3 and m q is the renormalization group invariant current quark mass. 
They are related to the scale dependent mass m q and the scale dependent chiral condensate (qq) A by the perturbative 
renormalization group equation 



(qq) = 



{qq)i 



[ln(A7A Q cD 2 )] a/2 ' 
m A [ln(A 2 /A QCD 2 )] a/2 



(3.4) 



(3.5) 



An overall factor (Nf = 2 times N c = 3) is omitted in Eq. (3.1). The chiral condensate (qq) and f v are known 
to be insensitive to the infrared regularization parameter p c [40]. Therefore we take p^/Aqcd 2 = e 01 and determine 
Aqcd to reproduce the pion decay constant = 93 MeV in the Pagels-Stokar formula [42] in the chiral limit. We 
obtain Aqcd = 738 MeV for Nf = 2 [10]. In the following calculations, we take A — 1 GeV in (3.5) and change the 
value of mJ =1GoV . For simplicity, we abbreviate m q =1GcY to m q below. 

At finite temperature and chemical potential, we use the imaginary time formalism [39] , and make the replacement 



(3.6) 



where ut n = (2n + l)7rT (n G Z) is the Matsubara-frequency for the quark. 10 

As a normalization, we define V by subtracting the c-independent part from V such that V reduces to the value 
of the free quark gas when a = 0. Sec Appendix A. 

We can study the chiral phase transition and the phase diagram by calculating V{<t 7 m q ] at given T and \x and 
searching the value of the chiral condensate ctq which minimizes the potential. The location of the first order phase 
transition line is determined by finding a gap in uq. In the chiral limit, do goes to zero smoothly as the second order 
phase transition line is approached from below. With finite quark masses, there is no distinct border between the 
symmetric and broken phases, and <io remains finite at all temperatures and chemical potentials. 

The phase diagram with several quark masses in the (T, /x) plane is shown in Fig. 3. The location of the tricritical 
point in the chiral limit is T t = 107 MeV and /it = 209 MeV. The open circles in Fig. 3 represent the critical end 
points for different quark masses. As shown in Fig. 4, the distance between TCP and CEP approximately scales as 
up to m q <~ 0(1) MeV, in agreement with (2.6). For larger masses, m q > 10 MeV, T c (m q ) does not change 
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much while n c (m q ) keeps on increasing. 
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FIG. 3. The phase diagram with several quark masses. The quark masses are evaluated at the momentum scale 1 GeV. 
The solid and dotted lines represent the first order and the second order phase transition, respectively. The filled circle is the 
tricritical point and open circles are the critical end points for different quark masses. 
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FIG. 4. The quark mass dependence of the critical temperature (upper figure) and the critical chemical potential (lower 
figure) . The slope of the solid line is | . 



B. The quark number susceptibility around CEP/TCP 

The quark number susceptibility \q is calculated from the normalized effective potential V as 

Fig. 5 and Fig. 6 are the results in the chiral limit. As we can see, \q 1S suppressed far below the chiral phase 
transition line and enhanced near TCP. In the chirally symmetric phase, Xq is equal to the value of the massless free 
quark gas x q Tee in this model. The region where Xq is enhanced is elongated in the direction parallel to the first order 
phase transition line. This is because the critical exponent for this direction (-y q = 1) is larger than for other directions 
("fq = |). We also found a jump in Xq along the second order phase transition line. Inside the critical region, however, 
the jump must be replaced by a cusp with certain critical exponents. See Appendix B. Our model can reproduce only 
the mean field behaviors. 
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FIG. 5. The quark number susceptibility near the tri-critical point in the chiral limit. The value of the susceptibility is 
divided by that of the massless free quark gas. The solid and dotted line represent the first and the second order phase 
transition, respectively and the filled circle is the tricritical point. 
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FIG. 6. The temperature dependence of \q a t fixed fj,'s. For /1 < 209 MeV, \ q has a jump across the second order phase 
transition line (0(4) line), which is consistent with the mean field theory. See Appendix B. 

Next we examine \ q for finite quark masses. Fig. 7 and Fig. 8. are the results for m q = 0.1 MeV and m q = 5.0 MeV, 
respectively. The location of CEP is (T c , /i c )=(104 MeV, 221 MeV) for m q = 0.1 MeV and (95 MeV, 279 MeV) for 
m q = 5.0 MeV. Xq diverges at CEP and is enhanced in the elongated region parallel to the first order phase transition 
line because the critical exponent is the largest for this direction as in the massless case. For m q = 0.1 MeV, TCP is 
still close to CEP and the elongated region includes the point (T t ,/i t ) while for m q — 5.0 MeV, the region deviates 
from it. 

At first sight, one might think that the analysis made in the previous section ceases to be valid at somewhere 
between m q = 0.1 MeV and m q = 5 MeV and the effect of TCP no longer survives for m q = 5 MeV which might be 



considered as the 'realistic' quark mass in this model. 11 However, this conclusion is too hasty. We will see in the 
next subsection that the hidden tricritical point still affects the physics near CEP even for m q = 5 MeV. 
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FIG. 7. The quark number susceptibility for m q = 0.1 MeV. The value of the susceptibility is divided by that of the massless 
free theory. The solid line is the first order transition line. The open circle represent the critical end point for m q = 0.1 MeV. 
The filled circle is at (T t ,/i t ). 
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FIG. 8. The quark number susceptibility for m q = 5 MeV. The value of the susceptibility is divided by that of the massless 
free theory. The open circle is the critical end point for m q — 5 MeV and the filled circle is at (T t , fit)- 



'In this model (qq) A=iGeV = ( — 276MeV) 3 at T = fi = 0. By using Gell-Mann-Oakes-Renner relation with m n — 140 MeV, 



m A=lGoV _ 4Me y^ 



C. The critical exponent for \q 



Now let us examine the critical exponent for Xq at CEP and TCP. We calculate it along the path parallel to the /j, 
axis in the T-fi plane from lower \i towards CEP/TCP at fixed T c or T t . 
First we consider the chiral limit. We expand V in the vicinity of TCP. 12 



V[a, m q = 0] = Vf rcc + 02 (T, [i)o 2 + a 4 (T, ^)a A + a 6 (T, M )cr 6 . 

The coefficients 02,04, 06 and Vf ree are summarized in Appendix A. cto is determined by the equation dV /da 
We obtain 



V[a = 0,TO 9 = 0] = Vf rco , 



above the chiral transition line, and 



V[a ,m g = 0] = Vfrcc + 7^2 ( 2a 4 _ 9a 2 a 6 ) - 7^3(04 - 3a 2 a 6 ) 3/2 , 



27aj 



27aj 



(3.8) 

= 0, 

TO 

(3.9) 
(3.10) 



below that line. \q i s obtained by taking the second derivative of (3.9) and (3.10) with respect to /1. Fig. 9 shows 
X q for numbers of |/i — yu t | 's. We determine the critical exponent j q defined in (2.14) numerically by using a linear 
logarithmic fitting 

lnx g = -7 9 m lM - Mtl + const. , (3.11) 
where const, is independent of [i. We obtained 7 9 = 0.51 ± 0.01 which is consistent with the mean field theory. 
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FIG. 9. The quark number susceptibility in the chiral limit as a function of \fj, — fj, t \ at fixed temperature T t . 



12 The reason for this expansion is twofold. First, in order to keep in line with the argument given in Section II. Second, 
technically we can approach TCP much closer to determine the exponent than directly reading it from Fig. 5. 



With finite quark masses, the expectation value cr is determined only numerically. This time we do not expand 
the potential around Co and directly read the exponent from Fig. 7 and Fig. 8. In Fig. 10, Xq is plotted for numbers 
of \n — /x c |'s for m q = 0.1, 5, and 100 MeV together with the calculated values of the critical exponent e defined in 
(2.9). 

ln% 9 = — eln \fi — fi c \ + const. . (3.12) 

For m q = 0.1 MeV we obtained e = 0.55 ± 0.02. This is significantly different from the prediction of the mean field 
theory e = |, which is a clear evidence of the effect of the tricritical point. We expect that the exponent changes 
towards | if we approach CEP much closer. 

For m q — 5 MeV, the slope of the data points changes at around \fi — fi c \ ~ 0.5 MeV. Therefore we fitted the data 
for \fi — fi c \ < 0.3 MeV and > 1 MeV separately and obtained the critical exponent 0.68 ± 0.02 for \fj, — fi c \ < 0.3 
MeV and 0.57 ± 0.01 for \fi — [i c \ > 1 MeV. We interpret this change of the exponent as the crossover of different 
universality classes discussed in the previous section. Note that the purely mean field-like exponent is seen in a very 
small region |/x — /i c \ < 1 MeV from CEP. This result is somewhat surprising to the present authors because, as seen 
in Fig. 8, TCP is far away from CEP already for m g = 5 MeV and the value of \q itself is unremarkable at (T t , fj, t ). It 
seems that, although the analysis in the previous section was made in the small quark mass limit, the effect of TCP 
is unexpectedly robust against the increase of the quark mass. 

As a check, we also calculated the exponent for m q = 100 MeV and obtained e = 0.64 ± 0.03 which is consistent 
with the mean field value |. For such a large quark mass, we see no indication of a change in the slope. The effect of 
TCP has completely disappeared. 
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FIG. 10. The quark number susceptibility for m q — 0.1, 5 and 100 MeV as a function of |p — u c \ at fixed temperature T c (m q ). 



IV. CONCLUSIONS 

Based on the universality argument and numerical model calculations, we studied the singular behavior of sus- 
ceptibilities near the critical/tricritical points. These two approaches are complementary, and we observed that the 
model calculation faithfully quantified the qualitative predictions obtained by the universality argument as long as 
the mean field behaviors are concerned. The important point is that, although we adopted a specific model, the 
qualitative behavior of \q is probably model independent. In particular, our results strongly suggest a possibility 
that the tricritical point affects the physics near the critical end-point. In other words, there are traces of the hidden 
tricritical point on the QCD phase diagram. Practically, the traces will be seen as the gradual change of the critical 
exponents since, after all, universality classes are characterized only by their critical exponents. It is expected that the 
exponents change from those of TCP to those of the Ising model via those of CEP in the mean field approximation. In 
order to really confirm this fascinating possibility, lattice simulations at finite chemical potentials [45] are necessary. 

Finally, we briefly comment on the implication of our results to heavy- ion experiments. The divergence of \q is 
directly related to an anomaly in the event-by-event fluctuation of baryon number B (divided by the entropy S) 

«^fffl (4.D 

which was originally introduced in [21] to probe the deconfined phase. Although neutrons are not observed, we expect 
that the event-by-event fluctuation of the proton number is relatively enhanced for collisions which have passed 'near' 
CEP/TCP. Pion and diphoton observables are discussed in [12-14]. As we remarked before, the critical exponents 
of the Ising model and the mean field theory are not so different numerically. Thus, the smallness of the critical 
region itself may not be an obstacle to the observability of critical phenomena in experiments. However, if we take 
the effect of TCP seriously either by assumption or stimulated by future lattice results, we must take into account the 
long- wavelength fluctuations of the pions as well as the sigma meson because the pions are no longer the 'environment' 
but participate in the critical fluctuations around the trace of TCP. 
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APPENDIX A: DESCRIPTION OF THE MODEL 

1. The normalized CJT effective potential 

We begin with the Cornwall- Jackiw-Tomboulis(CJT) effective potential [26] for QCD in the improved ladder ap- 
proximation [10] as a functional of the quark propagator S(p) at zero temperature and quark chemical potential after 
the Wick rotation, 

V[S] = V 1 [S] + V 2 [S}, (Al) 
Vi[S] = J ^tr {In [S^{p)S{p)] S^(p)S(p) + 1} , (A2) 



1 f f d 4 p d 4 k 



D^ip-k)}. (A3) 



Here "tr" is taken over the Dirac, flavor and color matrices (Gell-Mann matrices A°), and So(p) and D liv (p—k) are the 
free quark propagator and the gluon propagator in the Landau gauge (D tlv (p— k) = (5^ —Pfj,Pu/p 2 )/p 2 ), respectively. 
Vi[S] corresponds to the 1-loop potential with the quark 1-loop diagram and V^S 1 ] is the 2-loop potential with the 
one gluon exchange. 



We adopt the following approximation, the so-called Higashijima-Miransky approximation [37,38] for the QCD 
running coupling constant 



g 2 ((p - kf) - % 2 ^ V) + e{e P 2 )f(k") . 



(A4) 



where g is defined in (3.2). In this approximation with the Landau gauge, the renormalization of the quark wave 
function may be neglected at zero temperature and chemical potential. At finite temperature and chemical potential 
we need to take the wave function renormalization into account even in the Landau gauge [41]. However, we ignore 
this problem for the present purpose. Then the CJT effective potential can be rewritten as (3.1) in terms of the 
dynamical quark mass function £(p) using the corresponding Schwinger-Dyson equation for S(p). 

As a normalization, we define V by subtracting the cr-independent part from V such that V reduces to the value 
of the free quark gas when a = 0. We obtain 



V[(T,; 



Vfr, 



+ 



IT 2 ^ 

16T 

3C F air 2 



d\p\p In 



2 ,„ (S 2 (P 2 , u 2 n ; a, m q ) + p 2 + w 2 - M 2 ) 2 + 



(£ 2 (p 2 , ojI ■ 0, m q ) + P 2 +lo 2 - M 2 ) 2 + Vw 2 

p 2 + u>l)(p 2 + u>l + pi) [Hp 2 + lu 2 u + pD] 2 



+ pI) Mp 2 + ujI + p 2 c ) + p 2 + uj 2 



x m a a<j 



[ln(p 2 +^+pg)]~ 
(P 2 + ^1 +p 2 ) 3 

[Hp 2 



ln(p 2 +LJ 2 n +pl) + l 



■pDY 



ln(p 2 +^ 2 +p 2 ) + l-- 



(p 2 +C 2 +p 2 )4 

where a = 24/(llA c — 27V/) and the effective potential for the free quark Vfree is given by 

d 3 p 



21 



Vf r , 



-2T 



1 



(2tt) s 



(A5) 



(A6) 



In (l + e-^~^ T ) + In (l + e ~^' T ) 

with uj = ^Jp 2 +m 2 . In the chiral limit (m q = 0),the momentum integral can be easily performed and Vf ree becomes 

H 2 T 2 



Vfree {m q = 0) 



7tHT 4 



12tt 2 6 180 

The quark number susceptibility of the massless free quark gas is given by (omitting the overall factor NjN c ) 



x { r(m q = o) = 



(A7) 



(A8) 



2. The pion decay constant in the Pagels-Stokar formula 



The parameters p c and Aqcd are determined such that they reproduce the pion decay constant f„ = 93 MeV in 
the chiral limit. We calculate / w by the Pagels-Stokar formula [42] 



f 2 = 4A C 



d 4 p S(p;cr ,w g = 0) 
(2tt) 4 (E 2 (p;a ,m q = 0)+p 2 f 



Y,(p;a 0l m q = 0) 



p 2 a a ,m q = 0) 

2 dp 2 



(A9) 



In above equation, we set Nf — 2 because the pion consists of u and d quarks. 



3. Coefficients in the mean field expansion 

The explicit expressions of the coefficients 02, 04 and in Eq. (3.8) are 



(T 1 = 1tV (rl\ I J 2[ln(p 2 +^+^)] a - 2 (p 2 +^- M 2 ) 

a 2l i, M j ^^JWlP j (p2 + ^ +p 2 )2[(p2+LL) 2_ M 2 )2 + Vw 2 ] 

9 [ln(p 2 + + p 2 ) + 1 - f ] 2 (p 2 + ^ ) [ln(p 2 + + p 2 )] a ~ 2 
+ 2 (p 2 +^+p2)[( p 2 +c ,2 +?32)ln(p2+c , 2+p 2 )+p 2 +c ,2 ] > 

a (T ),i T V Adl I 2 [HP 2 + ^ + Pg)] 2 °' 4 [(p 2 + ^ - M 2 ) 2 - 

fT ) = _^ T ^ /di I 2 [MP 2 + ^ 2 +^)] 3a ' 6 (P 2 + ^-M 2 )[(P 2 4 

h' (p 2 +C 2 +p 2 )6[(p2 +LL) 2_ Ai 2 ) 2 



APPENDIX B: THE 0(4) CRITICAL LINE 

In this appendix, for completeness, we examine the singular behavior of \q along the 0(4) line emerging from TCP 
toward the temperature axis in the m = plane (See, Fig. 1). We call this line the 0(4) line because it consists of 
a sequence of critical points whose universality class is the same as that of the 0(4) spin model [2]. We again start 
with (2.1) with m = and the replacement a 2 -> 4> 2 = a 2 + (tt 1 ) 2 + (tt 2 ) 2 + (tt 3 ) 2 . The 0(4) line in the (T, fj,) plane 
is determined by the following equation 

a(2»=0. (Bl) 

Since b > does not vanish and smoothly varies along this line, we can drop the (j> 6 term. If we consider the mean 
field behavior, we can expand a around an arbitrary point (T c , fj, c ) on the line [31] 

a{T, l i) = C l {T-T c ) + D'{ l i- l i c ). (B2) 

In the mean field approximation, the (singular part of) thermodynamic potential becomes 

n MF = 0, (B3) 

above the 0(4) line, and 

a 2 

Vmf = (B4) 

below the 0(4) line. Taking the second derivative in ji, we see that the quark number susceptibility has a discontinuous 
jump across (T c , /j, c ) and that it is larger in the low temperature phase (below the 0(4) line) than in the high 
temperature phase (above the 0(4) line) except for points where D 1 = 0. Beyond the mean field approximation, we 
use the current theoretical estimate of the specific heat exponent of the 0(4) spin model [43] 

a 0.2. (B5) 

The minus sign means that the quark number susceptibility shows a cusp at T c as in the case of the A-point of liquid 
helium, [a is also negative for the 0(2) model.] Note that (T c , fj, = 0) is the point where D' = 0. It was shown in [44] 
that the 0(4) line is perpendicular to the temperature axis. Thus the quark number susceptibility has no singularity 
at (T c , l-i = 0) even in the chiral limit and increases monotonously as a function of the temperature, consistent with 
the results of lattice simulations. However, this smooth behavior is an exception only at fj, = 0. At any nonzero fi, 
X q has a cusp precisely at the critical temperature T c (/i). The cusp becomes higher and higher as we increase \i and 
finally diverges at the tricritical point. 
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